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SlMfcRY 


The  postbuckling  and  maximum  strength  analyses  of  uniformly  shortened, 
simply  supported  rectangular  plates  with  straight,  unloaded  edges  are  per- 
fcmed  by  using  Reissner' 6  variational  principle  in  conjunction  with  a 
deformation  theory  of  plasticity.  The  results  are  compared  to  (l)  ex¬ 
perimental  data  for  a  rectangular  plate  in  which  the  test  conditions 
basically  reflect  the  boundary  conditions  specified  in  the  present 
analyses,  and  (2)  a  potential  energy  solution  that  correlates  well  with 
experiment,  but  in  which  the  effect  of  waveform  change  on  the  average 
stress  carried  by  the  plate  is  accounted  for  only  in  a  gross  manner  and 
in  which  the  effect  of  small,  local  unloading  is  neglected.  Good  agree¬ 
ment  with  experiment  establishes  confidence  in  the  new  approach  and  in¬ 
dicates  that  the  simplified  technique  utilized  in  the  potential  energy 
solution  to  compensate  for  waveform  changes  may  be  employed  for  engineer¬ 
ing  purposes.  Both  analytical  approaches  predict  a  slightly  conservative 
plate  maximum  strength  relative  to  the  experimental  result.  This  dis¬ 
crepancy  is  attributed  to  slight  departures  of  the  test  conditions  from 
the  ideal  boundary  and  loading  conditions  assumed  in  the  analyses. 

An  application  of  Reissner's  principle  and  a  modified  version  of  the 
principle  is  undertaken  initially  to  obtain  the  elastic  postbuckling 
behavior  of  uniaxially  compressed  square  and  rectangular  plates.  Excel¬ 
lent  agreement  of  elastic-solution  results  with  essentially  "exact" 
solutions  of  other  authors  for  the  Bame  boundary  conditions  establishes 
the  effectiveness  of  the  Reissner  and  modified  Reissner  principles  and 
Justifies  the  application  of  the  Reissner-principle  approach  to  the 
maximum  strength  problem. 


iii 


TABLE  OF  CONTENTS 


SUMMARY  . 

LIST  OF  ILLUSTRATIONS  . 
LIST  OF  TABLES  .  .  .  . 
LIST  OF  SYMBOLS  .  .  .  . 


Page 

iii 

vi 

vii 

viii 


INTRODUCTION  . 

THEORY  . 

Statement  of  Problem  and  Basic  Assumptions  . 

Variational  Principle  . 

Specialization  for  Elastic  Behavior,  Unmodified  Version  .  .  . 
Specialization  for  Elastic  Behavior,  Modified  Version  .... 

METHOD  OF  SOLUTION  . 

Elastic  Problem,  Modified  Reissner  Functional  (Solutions  A)  .  . 
Elastic  Problem,  Unmodified  Reissner  Functional  (Solutions  B)  . 
Inelastic  Problem,  Unmodified  Reissner  Functional  . 

RESULTS  AND  DISCUSSION  . 

Elastic  Solutions  . 

Inelastic  Solutions  . 

CONCLUDING  REMARKS  . 

LITERATURE  CITED  . 

APPENDIXES 

I.  Euler  Equations  and  Boundary  Conditions  Derived  From  the 
Vanishing  of  the  First  Variation  of  the  Reissner  Func¬ 
tional  for  a  Two-Element  Plate  With  Prescribed  End 

Shortening  . 

II.  Modified  Reissner  Functional  for  Elastic  Problem  . 

III.  Unmodified  Reissner  Functional  for  Elastic  Problem  . 

IV.  Newton-Raphson  Iterative  Technique  . 

DISTRIBUTION  . 


1 

5 

5 

6 

9 

11 

13 

15 

17 

18 

22 

22 

23 


26 


40 


4? 

55 

59 

62 


6k 


LIST  OF  ILLUSTRATIONS 


Figure  ISSS. 

1  Plate  Jeometry  and  Coordinate  System  .  29 

2  Stress-Strain  Curve  for  Aluminum  2021-T3  Assumed  in 

Present  Inelastic  Analyses  .  30 

3  Load-Shortening  Curves  for  Square  Plates  Based  on 

Elastic  Solutions  .  31 

b  Load-Shortening  Curves  for  Rectangular  Plates  Based  on 

Elastic  Solutions  (b/X  =  1.25)  .  32 

5  Load-Shortening  Curves  for  Rectangular  Plates  Based  on 

Elastic  Solutions  (b/X  =  1.33)  .  33 

6  Load-Shortening  Curves  for  Rectangular  Plates  Based  on 

•  Elastic  Solutions  (b/X  =1.5)  .  3^ 

7  Load-Shortening  Curves  for  Rectangular  Plates  Based  on 

Elastic  Solutions  (b/X  =2.0)  35 

8  Comparison  of  Load-Shortening  Curves  for  a  Rectangular 
Plate  Based  on  Prediction  of  Present  Theory,  Analysis 

of  Reference  2,  and  Experimental  Results  of  Reference  3  •  36 

9  Theoretical  and  Empirical  Maximum  Strength  Criteria  for 

Simply  Supported  Rectangular  Plates  .  37 


vi 


LIST  OF  TABLES 


Table 

Page 

I 

Normalized  Stress  and  Displacement  Coefficients  for 
Square  Plates,  Solution  A1  . 

.  38 

II 

Normalized  Stress  and  Displacement  Coefficients  for 
Square  and  Rectangular  Plates,  Solutions  A2  . 

.  39 

III 

Normalized  Stress  and  Displacement  Coefficients  for 
Rectangular  Plates,  Solutions  A2  . 

40 

IV 

Normalized  Stress  and  Displacement  Coefficients  for 
Rectangular  Plates,  Solutions  B  (p  =  1.25  and  1.33)  .  • 

4l 

V 

Normalized  Stress  and  Displacement  Coefficients  for 
Rectangular  Plates,  Solutions  B  (p  =1.5  and  2.0)  ... 

.  42 

VI 

Normalized  Stress  and  Displacement  Coefficients  for 
a  Rectangular  Plate  (P  =  1.2  ,  ith/b  =  2.545  x  10-2)  . 
Inelastic  Results  Based  on  Aluminum  2024-T3  Plate 
Material . 

.  43 

VII 

Normalized  Stress  and  Displacement  Coefficients  for 
a  Rectangular  Plate  (P  =  1.  If  ,  Jth/b  =  2.545  x  10-2)  ; 
Inelastic  Results  Based  on  Aluminum  2024-T3  Plate 
Material  . 

.  44 

VIII 

Normalized  Stress  and  Displacement  Coefficients  for 
a  Rectangular  Plate  (p  =  1.6  ,  nh/b  =  2.545  X  10~2)  ; 
Inelastic  Results  Based  on  Aluminum  2024-T3  Plate 
Material . 

.  45 

vii 


LIST  OF  SYMBOLS 


A.  , B. .C 

1  1 


A.  . ,  B.  ..C.  . 

ij  ij  ij 

a.  . .b .  . . c.  . 
ij  10 


b 

D 


D 

s 


K 


E 


s 


e 


e 


cr 


F ' 


f 

h 


i.j.m.n 

K 


M  .M 
x-  y 

M 

xy 


m 


t 


in-plane  stress  parameters 
bending  stress  parameters 
in-plane  stress  parameters 
plate  dimension  in  y-direction,  in. 

Et3 

flexural  rigidity  of  solid  plate  =  - 5—  ,  lb-in. 

12(l-p  ) 

2Eh2tf  1 

flexural  rigidity  of  two-element  plate  =  -  ,  (u  =  g)  > 

lb-in.  3 

Young's  modulus  for  plate  material,  psi 

secant  modulus,  psi 

unit  shortening  applied  to  plate,  in. /in. 
buckling  strain  or  shortening,  in. /in. 

°eff 

complementary  energy  density  =  /  £effdtreff  psi 
displacement  parameter 

overall  thickness  of  two-element  plate  measured  between 
center  lines  of  faces,  in. 

integers 

material  constant  appearing  in  Equation  (7) 
numerical  constants 
plate  dimension  in  x-direction,  in. 
bending  moment  per  unit  length,  lb 
twisting  moment  per  unit  length,  lb 
number  of  buckles  in  x-direction 
materia],  constant  appearing  in  Equation  (7) 
thickness  of  solid  plate,  in. 

thickness  of  one  face  of  two-element  plate,  in. 


viii 


u 


11.  v.  w 


M  •  . 

1J 

V 


x.y.z 

a 

3 


xy 


xy 


xy 


eff 


e  .  c 

x  y 


✓  / 

€  ,  £ 

x’  y 


/  /  /  < 

£  .  € 

X  ■  y 


1 

X 


av 


cr 


cy 


eff 


Reissner  functional  .  lb-in. 


displacement  of  point  in  middle  surface  of  plate  in 
x-  .  y-  .  and  z-directions,  respectively,  in. 

nondimensionalized  displacement  parameters 

volume  of  plate,  in. 5 

plate  coordinates  (see  Figure  1  .  in. 

bending  stress  parameter  =  A^  ./A.  ^  -  31 3/^11  C13//Cll 

buckle  aspect  ratio  =  b/X 

total  shear  strain  in  xy-plane.  in. /in. 

shearing  strain  at  middle  surface,  in. /in. 

shearing  strain  due  to  twisting.,  in. /in. 


2 

effective  strain  =  - 


v/? 


in./ 


total  components  of  strain  in  x-  and  y-directions . 
respectively,,  in.  /in. 

components  of  strain  at  middle  surface  in  x-  and  y- 
directions.  respectively,  in. /in. 

components  of  strain  due  to  bending  in  x-  and  y- 
directions.  respectively,  in. /in. 

nondimensional  lateral  coordinate.  2y/b 

bucKle  half  wavelength  in  x-directicn.  in. 

Poisson's  ratio  for  plate  material 

nondimensional  axial  coordinate.  2x/X 

average  compressive  stress,  psi 

compressive  buckling  stress,  psi 

compressive  yield  stress  of  material,  psi 


effective  stress 


•v  2  2  Tl1 

\/<T  +  cr  -  cr  cr  +  it  .  psi 

V  x  y  x  y  xy  ^ 


■*  V 


tr 

max 


cr  ,cr 
x  y 


T 

XJr 

T  ' 

xy 
/  / 
‘xy 


maximum  average  compressive  stress,  psi 

components  of  stress  in  x-  and  y-directions, 
respectively,  psi 

local  average  stress  in  x-  and  y-direct^.ons, 
respectively,  psi 

bending  stress  in  x-  and  y-directions.  respectively,  psi 
total  shear  stress  in  xy-plane,  psi 
local  average  shear  stress,  psi 
shear  stress  due  to  twisting,  psi 


x 


INTRODUCTION 


The  analytical  deu^rmination  of  postbuckled  plate  behavior  in  the  plastic 
range  has  been  investigated  by  Mayers  and  Budiansky  for  square  plates 
in  Reference  1  and  by  Mayers,  Nelson,  and  Smith  for  rectangular  plates 
in  Reference  2.  The  latter  study  serves  to  explain  the  absence  of  a 
plate  maximum  strength  in  the  case  of  the  square-plate  analysis  of 
Reference  1  and  correlates  well  with  the  experimental  data  presented 
by  Stein  in  Reference  3.  In  Reference  4,  Stein  discusses  the  phenomenon 
of  buckle  wavelength  change  in  elastic  structures;  in  Reference  3,  a 
purely  theoretical  elastic  analysis  implies  that,  unlike  square  plates, 
postbuckled  rectangular  plates  can  change  buckle  aspect  ratio  at  reason¬ 
able  values  of  end  shortening.  The  experimental  evidence  presented  in 
Reference  3  proves  the  phenomenon  to  be  quite  true;  in  addition,  the 
experimental  data  reflect  a  maximum  load  which  cannot  be  predicted  by 
the  accompanying  elastic  analysis. 

Although  the  analysis  of  Reference  2  correlates  well  with  the  plate 
maximum  strength  data  of  Reference  3,  the  method  for  obtaining  the 
final  load-shortening  curves  leaves  a  question  unanswered.  This  ques¬ 
tion  concerns  the  validity  of  the  approximate  technique  of  accounting 
for  an  accurate  postbuckled  waveform  of  the  plate  by  adjusting  the 
first-approximation  load-shortening  curve  at  a  given  shortening  in  the 
plastic  range  by  the  ratio,  at  the  same  shortening  value,  of  the  "exact" 
elastic  load-shortening  curve  to  the  first -approximation  elastic  load¬ 
shortening  curve.  By  first-approximation  load-shortening  curve  is  meant 
that  curve  obtained  under  the  assumption  that  the  waveform  at  initial 
buckling  persists  through  the  postbuckling  range  to  plate  failure.  This 
technique  provides  good  correlation  of  the  theoretical  prediction  of 
Reference  2  for  the  maximum  strength  and  number  of  buckles  at  failure 
and  experimental  results  of  Reference  3.  This  procedure  is  highly 
desirable  for  engineering  purposes,  as  the  complexity  of  the  analysis 
is  reduced  considerably;  for  problems  in  which  an  accurate  elastic 
analysis  exists,  a  rapid  means  is  provided,  therefore,  for  evaluating 
the  maximum  strength  of  plate  elements.  Essentially  "exact”  elastic 
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solutions  to  the  problem  of  a  compressed,  simply  supported  flat  plate 
with  straight,  unloaded  edges  are  given  by  Stein  in  Reference  3  (square 
and  rectangular  plates)  and  by  Levy  in  Reference  5  (square  plates). 

One  way  in  which  to  justify  the  use  of  the  load -shortening  curve  adjust¬ 
ment  procedure,  on  other  than  an  experimental  correlation  basis,  is  to 
extend  the  potential  energy  solution  of  Reference  2  by  utilizing  suf¬ 
ficient  terms  in  the  functions  describing  the  displacements  u  ,  v  , 
and  v  such  that  the  essentially  exact  elastic  solution  can  be  obtained 
and  extended  into  the  inelastic  region.  This  procedure  would  not  require 
that  any  adjustment  be  made  to  the  inelastic  solution,  since  the  wave¬ 
form  beyond  initial  buckling  would  be  adequately  described. 

An  alternate  method  is  available,  however,  that  reduces  the  complexity 
of  the  problem  considerably.  This  method  is  based  on  the  variational 
principle  of  Reissner  given  in  Reference  6.  Stated  briefly,  the  prin¬ 
ciple  is  based  on  the  simultaneous  vanishing  of  the  first  variation  of 
the  Reissner  functional  with  respect  to  admissible  states  of  stress  and 
strain.  The  variation  yields  Euler  equations  corresponding  to  (l) 
stress-displacement  relations  and  (2)  equilibrium  conditions,  and  the 
associated  boundary  conditions.  Appendix  I  contains  the  variation  of 
the  Reissner  functional  for  a  two-element  plate  with  a  core,  rigid  in 
shear,  separating  the  faces.  This  configuration,  employed  in  References 
1  and  2  to  eliminate  the  added  complexity  of  plasticity  effects  through 
the  plate  thickness,  is  utilized  again  herein  for  comparison  purposes. 

Three  basic  advantages  of  employing  Reissner' s  principle  rather  than  the 
more  conventional  minimum  potential  energy  theorem  for  the  solution  of 
the  inelastic  plate  problem  are: 

1.  In  the  former,  stresses  and  displacements  are  assumed  as  com¬ 
pletely  independent  quantities  in  a  Rayleigh-Ritz  solution, 
while  in  the  latter,  displacements  must  be  determined  to  great 
accuracy  in  order  to  obtain  satisfactory  stresses  through  the 
particular  constitutive  relations  being  used. 

2.  A  first  approximation  to  the  inelastic  stress  distribution  can 
be  obtained  from  an  accurate  elastic  analysis;  thus,  the 


displacement  distributions  are  no  longer  required  to  high 
accuracy.  The  fact  that  the  stress -displacement  relations 
given  by  the  variation  of  the  Reissner  functional  are  satisfied 
only  approximately,  and,  perhaps,  unsatisfactorily,  by  a 
Rayleigh-Ritz  solution  is  of  little  consequence.  This  is  due 
to  the  proper,  independent  selection  of  stresses  and  displace¬ 
ments  which,  in  turn,  eliminates  the  need  to  evaluate  stresses 
in  terms  of  inaccurate  displacement  derivatives  through  the 
constitutive  relations.  With  equilibrium  and  compatibility  of 
deformations  satisfied  in  either  case  (potential  energy  method 
or  Reissner 's  principle),  it  would  appear  that  an  unacceptable 
stress  distribution  with  the  constitutive  relations  satisfied 
is  much  less  desirable  than  reasonably  correct  displacement  and 
stress  distributions  with  the  constitutive  law  not  satisfied 
exactly.  This  conjecture  is  strengthened  even  more  when  it  is 
considered  that  the  constitutive  law  assumed  to  hold  even  in 
elastic,  homogeneous,  isotropic  plate  theory  is  only  a  good 
approximation  for  a  conventional  material.  For  either  non¬ 
linear  elastic  or  inelastic  considerations,  the  desirability 
of  seeking  exact  satisfaction  of  an  idealized  constitutive  law 
should  be  subservient  to  insuring  that  displacement  compati¬ 
bility  and  equilibrium  are  satisfied. 

3.  The  effects  of  nonlinear  stress -strain  relations  are  readily 
incorporated  into  the  Reissner  functional  by  introducing  the 
Ramberg-Osgood  representation  of  the  material  stress-strain 
curve  (Reference  7)*  With  reference  to  the  remarks  in  item 
2,  it  should  be  noted  that  the  uniaxial  stress-strain  curve 
is  an  accurate  representation  in  general  only  up  to  the  region 
of  the  0.2$  offset  yield  stress. 

Strangely  enough,  Reissner 's  principle  has  not  been  widely  utilized  to 
date.  One  application  of  the  theorem  is  found  in  Reference  8,  in  which 
Sanders,  McComb,  and  Schlechte  extend  the  principle  to  consider  creep 
behavior  with  specific  applications  to  plates  and  columns.  In  Reference 
9,  Mayers  and  Rehfield  apply  a  modified  version  of  the  principle  to  study 
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the  postbuckling  behavior  of  axially  compressed  circular  cylindrical 
shells.  The  modification  of  the  principle  consists  of  coupling  the 
bending  stresses  or  their  resulting  moments  to  the  curvature  through 
Hooke's  law,  thus  leaving  only  the  in-plane  stresses,  in-plane  dis¬ 
placements,  and  lateral  displacements  as  variationally  independent 
quantities. 

The  present  analysis  also  employs  the  modified  principle  for  the  elastic 
postbuckled  plate  problem,  while  comparing  solutions  based  upon  two-  and 
three-term  lateral  displacement  functions  to  the  essentially  "exact" 
square-plate  solutions  of  Levy  and  Stein.  Both  the  modified  and  the 
unmodified  Reissner  functionals  are  employed  with  solutions  compared  to 
those  of  Stein  for  rectangular  plates  with  aspect  ratios  1.33,  1*5>  and. 
2.0  .  In  all  cases,  there  is  little  difference  obtained  in  the  quanti¬ 
tative  results  representing  the  load-shortening  relationships. 

With  confidence  established  in  the  present  method  of  solution,  the  in¬ 
elastic  postbuckling  solution  is  obtained  for  a  rectangular  plate  of 
aspect  ratio  5:1  .  Comparison  is  made  to  the  solution  of  Reference  2 
and  the  experimental  data  of  Reference  3.  Good  agreement  is  obtained, 
with  the  results  indicating  that  the  load-Bhortening  curve  adjustment 
procedure  employed  in  Reference  2  is  fully  Justified  for  predicting 
plate  maximum  strength  in  the  respective  analyses  of  References  1  and  2. 
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STATEMENT  CF  FROBIEM  AND  BASIC  ASSUMPTIONS 

The  problem  considered  herein  is  that  of  the  elastic  and  inelastic  post- 
buckling  behavior  of  uniformly  shortened,  simply  supported  square  and 
rectangular  plates  that  buckle  initially  in  the  elastic  range  and  that 
have  straight,  unloaded  edges  free  to  translate  in  the  plane  of  the 
plate.  As  observed  in  Reference  2,  major  factors  in  approximating 
satisfactorily  the  postbuckling  behavior  and,  ultimately,  the  maximum 
strength  of  uniformly  compressed,  simply  supported  rectangular  plates 
with  straight,  unloaded  edges  are  the  combined  effects  of  (l)  changes 
in  waveform  beyond  initial  buckling,  (2)  changes  in  buckle  aspect  ratio 
beyond  initial  buckling,  and  (3)  plasticity.  The  emphasis  is  placed  on 
determining  the  maximum  strength  of  rectangular  plates  without  utilizing 
the  load-shortening  curve  adjustment  procedure  of  References  1  and  2. 

The  adjustment  procedure,  in  gross  fashion,  accounts  for  change  in  buckle 
waveform  as  postbuckling  proceeds.  A  new  solution  of  the  maximum  strength 
problem  is  obtained  by  employing  Reissner's  variational  theorem;  the  re¬ 
sults  of  the  present  investigation  and  that  of  Reference  2  are  compared, 
and  the  accuracy  of  the  load-shortening  curve  adjustment  procedure  is 
evaluated.  To  establish  confidence  in  the  inelastic  analysis  based  upon 
the  Reissner  principle,  elastic  solutions  for  the  postbuckling  problem 
are  undertaken  and  compared  with  essentially  "exact”  solutions  presented 
in  References  3  and  5* 

The  inclusion  of  plasticity  effects  in  conjunction  with  analysis  based 
upon  the  small  strain,  moderate  rotations  strain-displacement  relations 
of  the  von  Karman  plate  theory  results  in  a  problem  that  is  nonlinear  in 
two  respects.  The  two-element  plate  description  (see  Figure  1),  introduced 
in  References  1  and  2  to  avoid  complexity  associated  with  accounting  for 
plasticity  effects  through  the  plate  thickness,  is  employed  herein  to 
provide  a  direct  comparison  of  the  results  of  the  present  analyses  with 
those  of  Reference  2.  The  core  separating  the  faces  of  the  plate  is 
considered  to  be  rigid  in  shear  in  the  present  investigation,  but  it 
provides  a  means  to  extend  the  current  analyses  to  include  transverse 
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shear  deformations  when  the  postbuckling  behavior  of  sandwich  plates  is 
of  interest.  Again,  plastic  behavior  is  based  upon  deformation  theory, 
and  the  ramifications  of  local  unloading  are  assessed  in  retrospect. 
Material  compressibility  is  neglected;  thus,  Poisson's  ratio  is  taken 
as  l/2  throughout  the  present  development. 

VARIATIONAL  FRINCIPIE 

In  essence,  Reissner's  principle  states  that  the  vanishing  of  the  first 
variation  of  the  Reissner  functional  with  respect  to  admissible  states 
of  stress  and  strain  establishes  the  stress-displacement  relations,  the 
equilibrium  equations  and  the  associated  boundary  conditions. 

The  general  functional  is  formulated  first  for  the  inelastic  problem; 
specializations  are  made  subsequently  for  purely  elastic  considerations 
in  both  an  unmodified  and  a  modified  version  of  the  functional.  The 
functional  modification  corresponds  to  the  one  introduced  in  Reference  9 
to  couple  the  bending  stresses  to  the  assumed  deflections  and  to  leave 
only  the  membrane  stresses  to  be  determined  from  the  variational  process. 

The  functional  for  plate -type  considerations  with  edge  displacements 
specified  is  given  as 
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a  e  + 
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where  F  '  is  the  complementary  energy  density.  If  the  end  loading  were 
to  be  prescribed  rather  than  the  end  shortening,  then  the  Reissner  func¬ 
tional  would  include  another  functional  term  corresponding  to  the  poten¬ 
tial  of  the  prescribed  edge  loads,  as  in  the  potential  energy  principle. 
The  vanishing  of  the  first  variation  of  the  functional  expressed  in 
Equation  (l)  is  performed  in  Appendix  I  for  the  two-element  plate  con¬ 
sidered  herein.  The  variation  establishes  Euler  equations  corresponding 
to  (l)  in-plane  and  out-of -plane  equilibrium  given  in  Equations  (35)  and 
(3  6),  respectively,  and  (2)  stress -displacement  and  moment -curvature 
relations  given  in  Equations  (34a)  and  (34b),  respectively.  Stress 
boundary  conditions  for  the  uniformly  compressed,  simply  supported  plate 
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with  straight,  unloaded  edges  resulting  from  the  variation  are  given  in 
Equation  (37  ). 

It  is  convenient  to  obtain  the  functional  of  Equation  (l)  in  terms  of 
stresses  and  strains  such  that  the  nonlinear  portion  of  the  functional 
due  to  plasticity  is  separable  from  the  elastic  portion.  To  this  end, 
the  complementary  energy  density  F '  is  first  written  incrementally  as 


dF '  =  e  do  +  e  do  +  7  dx 

tc  x  y  y  'xy  xy 


Next,  the  effective  stress  and  strain  quantities 
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which  are  related  by 
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are  introduced.  It  can  be  shown  readily  that 
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by  utilizing  the  secant  modulus  theory  stress-strain  relations  for  an 
incompressible  material  of  either  Reference  1  or  Reference  2.  The 
Reissner  functional  may  be  written  then  as 
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Equation '(6)  can  be  simplified  further  by  employing  the  Ramberg-Osgood 
three-parameter  representation  of  the  uniaxial  stress-6train  curve 
originally  developed  for  aluminum  alloy,  stainless  steel,  and  carbon- 
sheet  steel  in  Reference  7  and  given  by 


where  E  ,  K  ,  and  n  are  constants  determined  for  each  material  under 
consideration.  The  stress-strain  curve  (see  Figure  2)  for  aluminum 
202U-T3  used  in  the  present  inelastic  analysis  is  described  by 
K  »  3.1^  x  1017  ,  n  =  8.60  ,  and  E  =  10-7  x  10^  psi  .  Substitution 
of  Equation  (7)  into  Equation  (6)  and  subsequent  integration  yields 


Upon  integration  over  the  plate  thickness,  the  volume  integral  of 
Equation  (8)  for  the  two-element  plate  becomes,  with  V  =  2tfLb  , 


where  the  nondime  ns  ional  coordinates  ?  =  2x/X  and  tj  =  2y/b  are 
introduced,  and  the  total  stresses  and  strains  are  written  in  terms  of 
in-plane  and  bending  contributions  as  denoted  by  the  primed  and  double- 
primed  quantities,  respectively.  The  subscript  t  refers  to  the  top 
face,  and  the  subscript  b  refers  to  the  bottom  face.  Finally,  through 
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the  utilization  of  the  strain-displacement  relations  of  von  Karmen  plate 
theory,  modified  for  the  two-element  plate  as  indicated  in  Appendix  I, 
Equation  (27),  the  Reissner  functional  given  in  Equation  (9)  becomes 
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Specialization  for  Elastic  Behavior,  Unmodified  Version 

The  Reissner  functional  given  in  Equation  (l)  is  specialized  for  the 
elastic  postbuckling  problems  simply  by  setting  K  equal  to  0  .  The 

effective  stress,  in  general, 
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may  be  written  for  the  two-element  plate  as 
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It  is  easily  verified  that 
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Thus,  Equation  (lo),  with  K  equal  to  0  ,  can  he  written 
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The  formulation  of  the  Reissner  functional  given  in  Equation  (1I4-)  per¬ 
mits  an  independent  selection  of  the  stress  and  displacement  quantities 
that  satisfy  at  least  the  prescribed  displacement  boundary  conditions  im¬ 
plied  in  Appendix  I  for  application  of  the  Rayleigh-Ritz  procedure.  However, 
since  the  stress  distributions  are  selected  independently  in  the  Reissner 
formulation,  the  stress  boundary  conditions  may  be  satisfied  as  well. 

The  requirement  that  the  functional  be  stationary  with  respect  to  the 
free  parameters  which  describe  the  assumed  states  of  stress  and  displace¬ 
ment  then  approximately  satisfies  the  equilibrium  conditions  [Equations 
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(35)  and  (36)]  and  the  stress -displacement  and  moment -curvature  rela¬ 
tions  [Equation  (34)]. 

Specialization  for  Elastic  Behavior,  Modified  Version 

The  formulation  given  in  Equation  (14 )  is  modified  by  enforcing,  at  the 
outset,  the  moment- jurvature  relations  given  in  Equation  (34b),  which 
results  in  the  expressions 


Utilization  of  the  relations  given  in  Equation  (15)  permits  the  out-of- 
plane  equilibrium  Equation  (36 )  to  be  written 
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Equation  (16  )  is  analogous  to  the  familiar  fourth-order  out-of-plane 
equilibrium  equation  for  solid  plates  under  the  action  of  in-plane 
forces  given  by 
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Substitution  of  the  expressions  given  by  Equation  (15)  into  Equation  (14 ) 


[• 
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yields  the  modified  version  of  the  functional 


lr  -//|(r)[v+HM?)[^*H] 

+  +u  +w  w  l.ir/if +  /fatf  . 

\  E  /  L  ,X  'y  >X  >yl  2  L\e  /  Ve  / 

■  (rX,1) 

2 

♦  3f-*Tl+l[>  +v2  +W  V  +„2  1 

V  e  /  J  6  L  >**  »yy  >yy  >xy  J 

dldr)  (18) 

in  which  only  the  quantities  a'  ,  o'  ,  t'  ,  u  ,  v  ,  and  w  are  now 

x  y  xy 

variationally  independent. 
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METHOD  OF  SOLUTION 


The  most  efficient  application  of  Reissner's  variational  principle  to 
obtain  solutions  of  the  inelastic  postbuckled  plate  problem  utilizing 
a  direct  (Rayleigh-Ritz)  approach  requires  prior  knowledge  of  the  elas¬ 
tic  stress  and  displacement  distributions.  The  fact  that  the  plates 
considered  herein  buckle  elastically  suggests,  with  a  high  degree  of 
confidence,  the  use  of  stress  and  displacement  distributions  from 
existing  essentially  "exact"  elastic  plate  solutions  as  guides  in 
selecting  the  assumed  distributions  for  postbuckling  into  the  inelas¬ 
tic  range. 


In  general,  based  on  the  quite  accurate  elastic  stress  distributions 
depicted  in  Reference  5,  the  in-plane  stress  system  is  seen  to  be 
suitably  described  by  expressions  of  the  form 
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For  the  functional  given  in  Equation  (14 ),  the  bending  stresses  o"  , 
a"  ,  and  i''  are,  of  course,  expressible  in  a  similar  manner.  How- 

y  xy 

ever,  it  is  found  that  best  results  for  both  accuracy  and  simplicity 
are  obtained  in  the  present  problem  by  selecting  the  bending  stresses 
in  the  form  of  curvatures,  as  stipulated  by  Equation  (15 ),  but  with 
perfectly  free  amplitude  coefficients.  The  modified  functional  given 
in  Equation  (18  )  for  elastic  behavior  requires  no  specification  of  the 
bending  stresses,  since  the  moment-curvature  relations  are  enforced  at 
the  outset  through  Hooke's  law. 
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The  displacements  satisfying  all  geometric  boundary  conditions  (see 
Appendix  I)  are  suitably  expressed  by  series  of  the  form 

5  ■ <*•> 

even 

s  *  f©+n  <*»> 

even 

5  ■  1 1  v  «  ¥  -  *F  <**> 

odd 

where  e  is  the  applied  unit  shortening. 

The  coefficients  f  ,  u.  ,  v  ,  ,  a.  .  ,  b.  .  ,  and  c.  .  are  de- 

mn  mn  7  mn  ij  ij  ij 

termined  from  the  condition  that  5U"  =  0  with  respect  to  these  co¬ 
efficients.  The  number  of  free  variables,  indicated  by  Equations  (19) 
and  (20',  to  suitably  describe  the  stresses  and  displacements  may  seem 
excessive  at  first  glance;  however,  as  observed  in  Reference  9,  the  in¬ 
plane  stress  distribution  can  be  selected  such  that  one  or  both  of  the 
in-plane  equilibrium  equations  are  satisfied  independently  of  the  magni¬ 
tude  of  the  free  stress  coefficients.  If  the  y-direction  in-plane  equilib¬ 
rium  equation  (35b)  is  satisfied  independently  of  the  stress-coefficient 
magnitudes,  then  the  v-displacement  cannot  be  specified  arbitrarily.  This 
is  a  consequence  of  the  fact  that  the  satisfaction  of  the  y-direction  in¬ 
plane  equilibrium  equation,  which  is  performed  at  the  outset,  does  not 
permit  the  variation  B^U"  =  0  to  be  made.  Similarly,  the  u-displacement 
cannot  be  specified  in  terms  of  unknown  amplitudes  if  the  x-direction  in¬ 
plane  equilibrium  Equation  (35a)  is  satisfied  independently  of  the  stress- 
coefficient  magnitudes.  The  procedure,  wherein  one  or  both  of  the  in¬ 
plane  equilibriums  are  satisfied  independently  of  the  stress  coefficients, 
has  merit  in  problems  where  a  stress-function  approach  is  not  feasible. 

Such  is  the  case,  for  example,  in  Reference  9,  where  a  highly  nonlinear 
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elastic  shell  theory  is  employed,  and  in  the  present  problem,  where 
inelastic  behavior  is  present. 

ELASTIC  FROBIEM,  MODIFIED  KEISSNER  FUNCTIONAL  (SOLUTIONS  A) 


The  in-plane  stress  distribution  of  the  form  given  in  Equation  (19)  is 
assumed  initially,  for  the  present  analysis,  as 
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Due  to  the  large  number  of  stress  coefficients,  it  is  convenient  at  this 
point  to  invoke  the  stress  conditions  at  the  plate  edges  implicit  in 
Equation  (33)  and  thereby  reduce  the  number  of  free  variables  consid¬ 
erably.  The  satisfaction  of  the  stress  boundary  conditions  is  not 
actually  necessary,  since  the  surface  stresses  are  not  specified  for 
the  present  problem.  This  procedure  eliminates  8  stress  coefficients 
and  reduces  the  number  of  free  coefficients  to  10  .  In  view  of  the 
discussion  in  the  previous  section,  a  further  reduction  in  the  number 
of  free  stress  coefficients  can  be  made  by  taking  the  assumed  u-  and 
v-dis placement  functions  to  be  independent  of  free  parameters;  thus,  it 
is  possible  for  the  in-plane  equilibrium  Equation  (35)  to  be  satisfied 
Identically.  With  the  end  shortening  specified,  a  u-dis placement  func¬ 
tion  which  satisfies  the  geometric  boundary  conditions  and  which  contains 
no  free  parameters  is 
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Immediately,  Equation  (35*0  can  be  satisfied  identically  by  relating 

the  coefficients  of  o'  and  t' 

x  xy 

Now,  if  the  constant  v-dis placement  at  the  unloaded  edges  were  to  be 
specified  with  regard  to  magnitude,  then  the  function 

i  ■  '(f) 

would  satisfy  the  geometric  boundary  conditions  at  y  *  ±  b/2  and  would 

permit  the  coefficients  of  o'  and  t  '  to  be  related  such  that  Equa- 

y  xy 

tion  (35b)  could  be  satisfied  identically.  However,  although  the  v- 
dis placement  is  constant  on  the  unloaded  edges,  its  magnitude  f  is 
actually  unspecified  for  the  present  problem.  Indeed,  the  magnitude  of 
f  should  be  established  in  the  variational  process,  with  Equation  (35b) 
being  satisfied  only  approximately.  Nevertheless,  in  view  of  the  small 
influence  of  f  reflected  in  the  solutions  presented  in  References  1 
and  2,  it  is  found  that  the  specification  of  f  *  0  and  the  identical 
satisfaction  of  Equation  (35b),  in  addition  to  the  identical  satisfac¬ 
tion  of  Equation  (35a),  lead  to  accurate  results.  Thus,  Equation  (21) 
can  be  reduced  to  the  expressions  containing  only  six  free  variables 
given  by 


where  new  symbols  for  the  stress  coefficients  are  introduced  for 
simplicity. 


16 


There  are  no  restrictions  on  the  out-of -plane  displacement  w  other 
than  those  associated  with  the  geometric  boundary  conditions.  As  a  re¬ 
sult,  the  assumed  function  for  w  that  satisfies  the  geometric  boundary 
conditions  and  that  provides  sufficient  freedom  to  change  the  waveform 
as  postbuckling  progresses  is  obtained  by  retaining  only  the  w^  ,  w^1  , 
and  w^  terms  of  the  out-of -plane  displacement  series  given  in  Equation 
(20c).  The  actual  displacement  distributions,  with  both  in-plane  equi¬ 
librium  and  stress  boundary  conditions  satisfied,  are  assumed  throughout 
the  present  analyses  to  be 
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The  modified  Reissner  functional  given  in  Equation  (  18),  employing  the 
expressions  given  in  Equations  (22)  and  (23),  is  expanded  to  yield 
Equation  (39)  of  Appendix  II.  The  solution  obtained  by  utilization  of 
the  functional  given  in  Equation  (39)  with  (ith/b )  w^  =  w^  =  0  is 
designated  as  Solution  Al,  while  the  solutions  obtained  with  w^  re¬ 
tained  are  designated  as  Solutions  A2  . 


EIASTIC  iROBIEM,  UNMODIFIED  REISSNER  FUNCTIONAL  (SOLUTIONS  B) 

Instead  of  relating  the  bending  stresses  to  the  curvatures  through 
Hooke's  law,  as  is  the  case  in  Solutions  A,  the  bending  stresses  for 
the  unmodified  functional  expressed  by  Equation  (14)  are  selected  in¬ 
dependently.  With  the  utilization  of  the  three  out-of -plane  displace¬ 
ment  terms  given  in  Equation  (23),  the  bending  stress  distributions,  in 
the  form  of  the  curvatures  as  specified  by  Equation  (15),  satisfying 
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it  is  found  for  rec¬ 
tangular  plates  that  the  number  of  free  coefficients  may  be  reduced 
without  significant  loss  of  accuracy  by  setting  A^/A^  =  = 

C13/C11  =  a  and  A3i  =  B3i  =  C3i  =  0  *  Upon  substitution  of  the 

expressions  of  Equations  (22  ),  (23),  and  (24),  the  modified  Reissner 
functional  given  by  Equation  (14)  is  expanded  to  yield  Equation  (42)  of 
Appendix  III.  Solutions  obtained  on  the  basis  of  the  functional.  Equa¬ 
tion  (42),  with  w  given  by  Equation  (23)  are  designated  Solutions  B. 

INEIASTIC  PROBLEM,  UNMODIFIED  REISSNER  FUNCTIONAL 

For  the  determination  of  the  inelastic  postbuckling  behavior  and  maximum 
strength  of  plates,  the  Reissner  function  given  by  Equation  (10 )  is 
employed  in  conjunction  with  the  stress  distribution  described  by  Equa¬ 
tions  (22)  and  (24)  and  the  displacement  functions  of  Equation  (23). 
However,  due  to  the  Ramberg-Osgood  generalized  stress-strain  curve 
parameter  n  being  large  for  conventional  metals  in  general,  the  in¬ 
dicated  integration  in  closed  form  of  the  effective  stress  terms  in 
Equation  (10)  is  prohibitive.  Consequently,  the  integration  is  per¬ 
formed  numerically. 


the  boundary  conditions  [Equation  (37)],  are 
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where  all  amplitude  parameters  are  free;  however, 
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For  K  =  0  ,  Equation  (10)  reduces  to  the  functional  given  by  Equation 

(42).  Hence,  the  functional  for  the  inelastic  problem  may  be  written  as 
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where  (U' //®V)eiastic  S*ven  in  Appendix  III  by  Equation  (4 2) . 

The  required  solutions  are  obtained  by  determining  the  values  of  the 
stress  and  displacement  coefficients  that  render  the  formulated  Reissner 
functionals  stationary  in  value.  To  achieve  elastic  Solutions  Al,  A2, 
and  B  and  the  inelastic  solution,  the  procedure  requires  the  solutions 
of  8,  9>  13,  and  13  simultaneous  nonlinear  algebraic  equations  re¬ 
spectively.  The  method  utilized  in  the  solution  of  these  equations  is 
a  basic  Newton-Raphson  iteration  procedure.  Appendix  IV  contains  a  de¬ 
tailed  description  of  the  method.  The  modified  steepest  descents  min¬ 
imization  technique  developed  in  Reference  1  and  improved  in  Reference  2 
to  obtain  relative  minima  of  the  strain-energy  functional  is  not  as 
efficient  for  the  present  solutions,  since  the  vanishing  of  the  first 
variation  of  the  Reissner  or  modified-Reissner  functionals  with  respect 
to  purely  arbitrary  stress  and  displacement  states  does  not  correspond 
to  the  achievement  of  relative  minima. 


The  application  of  the  Newton-Raphson  technique  yields  numerical  values 
for  the  coefficients  which  can  be  utilized  in  conjunction  with  the  axial 
stress  distribution  to  construct  load-shortening  curves.  The  load¬ 
shortening  curves  are  plots  of  average  axial  stress  versus  unit  shorten¬ 
ing,  normalized  with  respect  to  the  buckling  stress  and  strain  for  a 
square  (or  infinitely  long),  simply  supported  plate,  respectively.  The 
average  stresses  in  the  present  analyses  are  obtained  directly  from  an 
integration  of  a ^  over  the  plate  width.  Numerical  values  of  the  stress 
and  displacement  coefficients,  normalized  with  respect  to  the  critical 
strain  of  a  simply  supported  square  (or  infinitely  long  plate),  are 
listed  in  Tables  I  through  VIII  for  various  values  of  unit  shortening. 
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Two  solutions,  denoted  by  Solution  A1  and  Solution  A2,  are  effected  for 
square  plates  employing  the  modified  elastic  functional  and  are  compared 
in  Figure  3  to  the  first-approximation  elastic  solution  of  Reference  2 
and  the  essentially  "exact"  solutions  of  References  3  and  5;  corresponding 
normalized  stress  and  displacement  coefficients  for  specific  values  of 
unit  shortening  are  listed  in  Tables  I  and  II.  Solution  Al,  which 
utilizes  only  the  and  w^  terms  in  the  wniisplacement  function, 

is  accurate  only  for  square  plates.  Analyses  results  for  rectangular 
plates  with  aspect  ratios  1.25  ,  1.33  ,  1.5  >  and  2.0  ,  obtained  from 
utilization  of  the  modified  (Solutions  A2)  and  unmodified  (Solutions  B) 
elastic  functionals,  are  compared  in  Figures  4  througn  7  to  the  first- 
approximation  elastic  solutions  of  Reference  2  and  the  essentially  "exact" 
solutions  of  Reference  3.  Tables  II  through  V  contain  representative 
values  of  the  normalized  stress  and  displacement  coefficients,  at  given 
unit  shortenings,  for  these  cases. 

A  maximum  strength  analysis  of  an  aluminum  2024-T3  rectangular  plate 

with  aspect  ratio  5  ,  for  which  test  data  exist  in  Reference  3,  is 

obtained  through  utilization  of  the  unmodified  inelastic  functional. 

The  two-element  plate  geometry  analyzed  is  made  effectively  equivalent 

to  the  plate  tested  in  Reference  3  by  equating  the  expression  for  the 

critical  stress  of  the  two-element  plate,  derived  in  Reference  2  to  be 

o  *  4/3[E(nh/b)  ]  ,  to  the  classical  critical  stress  expression 

2  2 

ocr  *  E(irt/b)  /3(l-p  )  for  a  solid  plate.  The  load-shortening  curve 
obtained  from  this  analysis  is  compared  in  Figure  8  to  the  test  data 
of  Reference  3  and  the  analysis  results  of  Reference  2.  Tables  VI 
through  VIII  contain  numerical  values  for  the  elastic  and  inelastic 
stress  and  displacement  solution  coefficients,  at  various  values  of 
unit  shortening,  corresponding  to  the  buckle  wavelengths  considered. 

The  determination  of  the  values  of  unit-shortening  ratios  at  which  the 
buckle  wavelengths  change  as  the  unit  shortening  increases  is  described 
in  detail  in  Reference  2  and  is  briefly  restated  herein  for  completeness. 
The  calculation  is  made  by  construction  of  load- shortening  curves  cor¬ 
responding  to  each  buckle  half  wavelength  X  or  the  number  of  buckles  m  . 
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For  a  given  plate  geometry,  the  load-shortening  curves  corresponding  to 
m  and  nri-1  buckles  intersect.  At  a  particular  value  of  shortening, 
the  Reissner  functional  magnitude  (area  under  the  load-shortening  curve) 
corresponding  to  mtl  buckles  will  be  equal  to  that  corresponding  to 
m  buckles.  Theoretically,  this  value  of  unit  shortening  is  the  point 
at  which  a  "Jump"  is  first  possible.  For  additional  shortening,  the 
curve  corresponding  to  bh-1  buckles  becomes  temporarily  the  appropriate 
load-shortening  path.  The  procedure  is  followed  with  successively  in¬ 
creased  numbers  of  buckles  considered  to  the  point  at  which  a  maximum 
load  is  reached. 


21 


RESULTS  AMD  DISCUSSION 


EIASTIC  SOLUTIONS 

To  establish  confidence  in  the  method  of  analysis  enployed  herein,  ap¬ 
plications  of  Reissner's  principle  and  a  modified  version  of  the  prin¬ 
ciple  have  been  made  for  comparison  with  essentially  "exact"  solutions 
to  the  elastic  postbuckling  problem  of  uniaxially  compressed,  simply 
supported  plates  possessing  straight,  unloaded  edges.  Excellent  agree¬ 
ment  of  the  present  results  with  the  essentially  "exact"  solutions  of 
Stein  (Reference  3)  and  Levy  (Reference  5),  shown  in  Figures  3  through 
7,  indicates  that  the  method  of  solution  is  sufficiently  flexible  to 
permit  appreciable  simplifications  to  be  made  without  significant  loss 
of  accuracy. 

One  simplification,  that  of  coupling  the  bending  stresses  to  the  curva¬ 
tures,  introduced  in  Reference  9  and  employed  herein  to  perform  Solu¬ 
tions  A,  is  seen  to  be  a  desirable  procedure  for  elastic  plate-type 
problems,  in  that  the  number  of  free  stress  parameters  in  the  direct 
variational  solution  is  considerably  reduced.  Another  simplification, 
that  of  describing  the  u-displacement  function  by  a  single  term  and  the 
specification  of  v  =  0  (utilized  in  both  Solutions  A  and  B),  which 
allows  the  in-plane  equilibrium  equations  to  be  satisfied  identically, 
appears  to  be  justified  by  the  excellent  agreement  obtained  upon  com¬ 
parison  with  the  essentially  "exact"  solutions.  Obviously,  if  the 
stresses  were  required  to  be  found  in  terms  of  the  displacements  through 
application  of  the  constitutive  law  (potential  energy  approach),  the  in¬ 
fluence  of  the  simplified  u-  and  v-dis placement  functions  would  have  a 
profound  effect,  relative  to  the  satisfaction  of  in-plane  equilibrium. 
Hence,  the  feature  of  selecting  the  stresses  and  displacements  inde¬ 
pendently  and  of  satisfying  Hooke's  law  only  approximately,  as  per¬ 
formed  herein  for  Solutions  B  with  the  aid  of  Reissner's  variational 
principle,  offers  a  very  desirable  approach  relative  to  that  of  po¬ 
tential  energy,  wherein  very  accurate  displacement  distributions  are 
required  to  attain  satisfactory  stress  distributions  through  the  con¬ 
stitutive  law. 
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The  reduced  complexity  of  the  problem,  as  a  consequence  of  the  invoked 
simplifications,  is  strikingly  illustrated  by  the  fact  that  Levy 
(Reference  5)  requires  at  least  four  out-of -plane  displacements  and 
fourteen  stress  function  terms  to  obtain  the  essentially  "exact" 
solution  shown  in  Figure  3  for  the  square  plate  at  the  maximum  shorten¬ 
ing  ratio  considered.  At  the  shortening  ratio  corresponding  to  the 
calculation  limit  of  Stein's  solution  (Reference  3),  at  least  three 
out-of -plane  displacement  terms  and  either  ten  stress -function  terms 
or  ten  in-plane-displacement  terms  are  required  to  effect  the  "exact" 
solution  by  the  Levy  and  Stein  approaches,  respectively. 

IMEIASTIC  SOLUTION 

The  maximum  strength  analysis  of  a  uniformly  shortened,  simply  supported 
rectangular  plate  based  on  Reissner's  variational  principle  has  been 
performed  to  Justify  the  procedure  for  adjusting,  in  a  gross  manner, 
the  load-shortening  curve  to  account  for  waveform  changes  beyond  ini¬ 
tial  buckling.  This  procedure  was  utilized  in  the  solutions  of 
References  1  and  2  in  conjunction  with  the  potential  energy  method  and 
with  the  assumed  persistence  during  postbuckling  of  the  waveform  oc¬ 
curring  at  initial  buckling.  The  load -shortening  curve  obtained  in  the 
present  analysis  for  the  particular  plate  geometry  reflected  in  the 
experimental  work  of  Reference  3  agrees  well  with  the  results  of 
Reference  2,  as  shown  in  Figure  8.  However,  both  the  present  results 
and  the  results  of  Reference  2  are  slightly  conservative  with  respect 
to  the  experimental  data  of  Reference  3;  the  increased  slope  of  the 
experimental  curve  at  initial  buckling  indicates  that  the  test  data 
have  been  influenced  to  some  degree  by  the  test  fixtures. 

To  place  the  present  maximum  strength  solution  in  more  general  perspec¬ 
tive,  it  is  plotted  in  Figure  9  along  with  a  portion  of  the  theoretical 
and  experimental  plate  maximum  strength  data  compiled  in  Figure  5  of 
Reference  2.  The  reference  numbers  in  Figure  9  are  those  of  the 
present  report.  It  can  be  seen  that  the  maximum  strength  analyses  of 
both  References  2  and  3  correlate  more  closely  with  the  test  data  of 
Reference  11  than  with  the  test  point  taken  from  Reference  3.  However, 


as  mentioned  previously,  the  test  data  of  Reference  3  reflect  an  appar¬ 
ent  increase  in  plate  strength  due  to  test  fixture  effects;  however, 
since  the  actual  cause  of  the  increased  slope  of  the  load- shortening 
curve  at  initial  buckling  (see  Figure  8)  could  not  be  isolated,  no 
attempt  has  been  made  to  adjust  the  experimental  postbuckling  and  maxi¬ 
mum  strength  results.  The  comparison  of  the  present  results  and  those 
of  Reference  2  with  the  maximum  strength  criterion  based  upon  the 
von  Karman  effective -width  concept  (Reference  10 )  suggests  that  when 
the  buckling  stress  for  a  simply  supported  plate  is  greater  than 

0.3  o  ,  the  curve  faired  through  the  three  points  given  by  the 

cy 

calculated  point  of  the  present  analysis  and  the  calculated  points  of 

Reference  2  at  (acr)c]/°cy  =  0*3  and  0.6  ,  respectively,  should  be 

utilized  for  maximum  strength  prediction  of  rectangular  plates 

(L/b  >5)  •  Application  of  the  present  analysis  technique  to  improve 

the  results  of  Reference  2  at  (a  )  ,/a  =  0.3  and  0.6  ,  respec- 

tively,  would  show  little  effect,  since  the  higher  initial  buckling 

stresses  reduce  the  postbuckling  range  considerably  prior  to  the 

occurrence  of  a  maximum  load.  In  fact,  a  glance  at  Figure  5  of 

Reference  2  shows  that  the  maximum  strengths  of  simply  supported  plates 

with  various  types  of  unloaded-edge  conditions  after  buckling  reflect 

little  difference  as  (a  )  ,/a  approaches  about  0.6  . 

v  cr,cl/  cy 

The  complexity  of  extending  the  potential  energy  (Rayleigh-Ritz) 
solution  of  Reference  2  (without  the  load-shortening  curve  adjustment 
procedure)  to  obtain  results  equivalent  to  those  of  the  present 
analysis  is  prohibitive,  due  to  the  number  of  displacement  terms  re¬ 
quired  to  obtain  an  accurate  in-plane  stress  distribution  through  the 
constitutive  law.  The  procedure  employed  in  References  1  and  2  re¬ 
quires  the  existence  of  an  "exact"  elastic  solution  upon  which  to  base 
the  adjustment  procedure;  on  the  other  hand,  efficient  use  of  the 
present  theory  requires,  in  the  absence  of  stress  distribution  infor¬ 
mation  determined  on  the  basis  of  limited  experimental  data,  no  more 
than  approximate  elastic  solutions  from  which  forms  of  the  assumed 
stress  distributions  can  be  obtained. 
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Another  advantage  of  the  present  method  of  analysis  relative  to  the 
potential  energy  approach  of  References  1  and  2  is  that  in  the  former, 
the  stress  distribution  can  be  evaluated  throughout  the  plate  after  the 
stationary  value  of  the  Reissner  functional  has  been  determined,  whereas 
in  the  latter,  the  stress  distribution  is  known  only  for  the  solution 
based  upon  the  initial  buckling  waveform.  The  stress  distributions 
corresponding  to  the  adjusted  load-shortening  curves  remain,  of  course, 
unknown. 

Due  to  the  utilization  of  a  deformation  theory  of  plasticity,  occurrence 
of  unloading  in  particular  regions  of  the  plate  as  the  shortening  and 
number  of  buckles  both  increase  must  be  considered.  It  is  found  in  the 
present  analysis  that  local  unloading  of  small  magnitude  (no  greater 
than  5%)  does  occur  in  inelastic  regions  of  the  plate  when  the  buckle 
wavelength  is  assumed  to  change  abruptly;  however,  due  to  the  known 
inaccuracy  of  the  uniaxial  stress -strain  curve  representation  for  a 
material  and  taking  cognizance  of  the  fact  that  the  total  load  carried 
by  the  plate  is  an  Integrated  phenomenon,  the  small,  local  unloading 
that  does  occur  in  the  inelastic  portion  of  the  plate  is  ignored  in  the 
present  problem.  Thus,  it  may  be  considered  that  the  present  problem 
has  been  undertaken  under  the  assumption  that  the  material  is  simply 
of  a  nonlinear  elastic  nature. 

The  simplifying  assumption  in  the  present  solution  that  the  v -displacement 
is  zero  throughout  the  plate  should  be  assessed.  Based  on  the  correlation 
of  the  results  for  the  elastic  solutions  carried  out  herein  and  in  Ref¬ 
erences  3  and  5,  the  influence  of  the  violation  of  the  condition  that 
the  unloaded  edges  be  free  to  translate  is  considered  to  be  minor.  The 
actual  correlation  is  shown  in  Figures  3  through  7. 
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CONCLUDING  REMARKS 


The  basic  factors  required  to  accurately  predict  the  Inelastic  post- 
buckling  behavior  and  the  maximum  strength  of  simply  supported  rec¬ 
tangular  plates  with  straight,  unloaded  edges  have  been  shown  to  be 
(l)  buckle  wavelength  changes,  (2)  changes  in  waveform  beyond  initial 
buckling,  and  (3)  the  onset  of  plasticity.  The  impact  of  item  (l)  can 
be  seen  in  the  inelastic  postbuckling  analysis  of  simply  supported 
square  plates  presented  in  Reference  1.  While  the  buckle  aspect  ratio 
remains  unity,  no  clearly  defined  plate  maximum  strength  is  obtained 
for  the  quite  reasonable  range  of  plate  end  shortening  considered. 
However,  the  consideration  of  all  three  factors  collectively  for  the 
first  time  in  the  rectangular  plate  investigation  of  Reference  2  re¬ 
sults  in  an  analytical  solution  for  the  prediction  of  the  maximum 
strength  of  a  simply  supported  plate,  which  correlates  excellently 
with  experimental  data  presented  in  Reference  3*  The  experiment  re¬ 
flects  loading  and  boundary  conditions  very  close  to  those  assumed  in 
analyses  of  References  1  and  2.  With  regard  to  item  (2),  the  analyses 
of  References  1  and  2  utilize  a  load-shortening  curve  adjustment  pro¬ 
cedure  that  only  grossly  accounts  for  the  waveform  changes  beyond  in¬ 
itial  buckling  as  they  effect  the  load-shortening  curve.  The  present 
analysis,  based  upon  Reissner's  variational  principle,  establishes  the 
validity  of  the  adjustment  procedure  and,  in  addition,  illustrates 
another  attractive  method  for  obtaining  engineering  solutions  to  plate 
maximum  strength  problems.  As  concerns  item  (3),  the  ability  of  the 
secant  modulus  deformation  theory  of  plasticity  to  adequately  represent 
inelastic  behavior  is  demonstrated.  Local  unloading  of  small  magnitudes 
occurs  in  certain  inelastic  portions  of  the  plates  analyzed;  however, 
the  influence  of  small  unloading  on  the  postbuckling  and  maximum  plate 
strengths  (which  are  integrated  phenomena)  is  not  significant. 

Efficient  use  of  the  potential  energy  method  of  solution  with  the 
adjustment  procedure  incorporated  requires  that  essentially  exact 
elastic  solutions  be  available  to  provide  the  necessary  basis  for 
carrying  out  the  adjustment  procedure.  However,  should  neither  "exact ' 


26 


elastic  solutions  exist  nor  limited  experimental  data  be  available  to 
suggest  the  form  of  the  stress  distributions,  then  only  approximate 
elastic  solutions  are  required  to  obtain  the  general  form  of  the  dis¬ 
tributions  to  be  used  in  conjunction  with  Reissner's  variational 
principle. 

The  fact  that  Reissner's  principle  has  had  little  application  to  date 
may  be  attributed  to  the  apparent  complexity  involved  in  the  require¬ 
ment  to  assume  the  states  of  stress  and  displacement  independently. 
However,  the  present  analyses  indicate  that  the  principle  has  suf¬ 
ficient  flexibility  to  permit  major  simplifications  to  be  made  with 
little  loss  in  accuracy.  In  particular,  the  features  of  (1)  employing 
only  the  first  approximation  of  the  in-plane  displacements  and  (2) 
coupling  of  the  bending  stresses  to  the  curvatures  through  the  con¬ 
stitutive  law  in  the  elastic  case  (the  modified  Reissner  principle 
introduced  in  Reference  9)  are  found  to  be  extremely  useful  in  re¬ 
ducing  the  number  of  free  amplitude  parameters  appearing  in  the  assumed 
stress  and  displacement  functions.  Item  (1),  although  resulting  in  a 
violation  of  the  unloaded -edges  boundary  condition  and  permitting  only 
approximate  satisfaction  of  the  constitutive  law,  is  seen  to  have  only 
minor  influence  on  the  accuracy  of  the  elastic  solutions.  The  infer¬ 
ence  of  item  (2)  is  also  found  effective  for  the  present  inelastic 
solution,  where  the  bending  stresses  are  assumed  in  the  form  of  the 
curvatures  as  indicated  by  the  constitutive  law  but  with  completely 
free  amplitude  coefficients  whose  magnitudes  are  determined  from  the 
variational  process. 

The  unconservative  nature  of  the  von  Karmen  effective-width  formula, 

when  applied  to  simply  supported  plates  possessing  a  buckling  stress 

greater  than  0.3  (o  )  ,  suggests  that  the  combined  results  of  the 

cy 

present  inelastic  analysis  and  the  inelastic  analyses  of  Reference  2 

be  utilized  for  design  purposes.  However,  additional  tests,  utilizing 

the  technique  described  in  Reference  3,  are  required  for  the  range 

0.3  <  (°cr)ci/°cy  <  0»6  .  For  simply  supported  plates  with  buckling 

stresses  less  than  0.3  (o  )  ,  the  von  Kdrwdn  effective-width  formula 

cy 


2" 


provides  conservative  design  data  for  plates  in  which  the  edges  do  not 
move  out  of  the  original  plane  in  the  postbuckling  range.  Where  other 
plate  loading  and  boundary  conditions  are  of  interest,  the  techniques 
of  either  Reference  2  or  the  present  inelastic  analysis  provide  an 
efficient  means  for  obtaining  maximum  strength  behavior.  The  com¬ 
bination  of  the  results  of  such  analyses  with  a  minimum  of  Judicious 
testing  (for  correlation  purposes)  could  provide  a  useful  tool  for 
removing  some  of  the  conservatism  from  present  design  analysis 
techniques . 
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STRESS,  KSI  (TENSION  OR  COMPRESSION 
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Figure  2.  Stress-Strain  Curve  for  Aluminum  202U-T3  Assumed  in  Present 
Inelastic  Analyses. 
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Figure  L.  Load-Shortening  Curves  for  Rectangular  Plates  Based  on  Elastic 
Solutions  (b/X  =  1.25)  . 
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Figure  5.  Load-Shortening  Curves  for  Rectangular  Plates  Based  on  Elastic 
Solutions  (b/X  =1.33)  . 
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Figure  6.  Load-Shortening  Curves  for  Rectangular  Plates  Based  on  Elastic 
Solutions  {b/\  =1.2)  . 
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Figure  7-  Load- Shortening  Curves  for  Rectangular  Plates  Based  on  Elastic 
Solutions  (b/X  =  2.0)  . 
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NUMBER  OF  BUCKLES 
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Figure  8.  Comparison  of  Load-Shortening  Curves  for  a  Rectangular  Plate  Based  on  Prediction  of 
Present  Theory,  Analysis  of  Reference  2,  and  Experimental  Results  of  Reference  3. 
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NORMALIZED  STRESS  AND  DISPIACEMEKT  CCEFFICIEHTS 
FOR  SQUARE  PIATES,  SOLUTION  A1 
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Ai'PjvJDlX  : 

EULER  EQUATIONS  AND  BOUNDARY  CONDITIONS  DERIVED  FROM  THE  VANISHING 
OF  THE  FIRST  VARIATION  CF  THE  RE  ISSUER  FUNCTIONAL  FOR 
A  TWO-EIEME17J  PIATE  WITH  PRESCRIBED  END  SHORTENING 

The  Reissner  functional  for  present  considerations  is  defined  as 


U"  =  /  lag  +  a  e  +  t  7  -F'l  dV  (?6) 

J  1  x*x  y  y  xy'xy  |  v 

V 

The  strain-displacement  relations  for  small-strain,  moderate  rotations 
plate  theory,  when  modified  for  the  two-element  plate  of  Figure  1,  be¬ 
come 


1  2 
o  v 

2  ,x 


1  2 
o  v 

2  ,y 


w 

,yy 


)  (27) 


y'  ±  y"  = 
xy  'xy 


w  w  ±  hw 
,x  ,y  ,xy 


where  the  +  and  -  signs  correspond  to  the  top  and  bottom  elements 
of  the  plate,  respectively. 


The  stresses  may  be  written 


a'  ±  a" 

y  y 


T  " 

xy 


±  T 


xy 


)  (28) 


J 


where  the  primed  and  double-primed  quantities  denote  in-plane  and 
bending  contributions,  respectively. 
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Substitution  of  Equations  (27)  and  (28)  into  Equation  (26)  yields 


The  average  stresses  and  resulting  bending  moments  are 


hQ 


I 


(Continued) 


M  =  t  -  £  T  -T  =  t-hT  '  ' 

xy  f  2  ^  xyt  xyb  I  f  xy 


Substitution  of  Equation  (30)  into  Equation  (20)  and  with  (F^  +  F^) 
replaced  by  2F'  gives 


U"  =  t. 


L/2  b/2 


-L/2  -b/2 


2o '  (  u  +  i  w2  ^  +  2c  '  (  v  +iv2\ 
x  ^  ,x  2  ,x  J  y\  ,y  2  ,y  J 


+  2t  '  u  +v  +ww  ] - w 

xy  ^  ,y  ,X  ,x  ,y  j  ,xx 

M  M 

-  -2-  w  +2  w  -  2F  '  dxdy 

t„  *yy  >xy 


In  the  elastic  case, 


1  2  2 
-  2E  |%  +  V 


1  2  2  2  2 
-acr  +3t  +a  +  a  -  a  a  +  3* 

t  xt  xyt  *b  yb  *b 


■'.ip 

2E  LX 


2  +  a  ^  -  a  'a '  +  3t  '2  +  - =•  (/  +  1^  -  M  M  + 

y  X  y  Xy  (tfh)2  X  y  X  y 


3^  )1 

xy  j 


7eff 

whereas  in  the  inelastic  case,  F'  =  /  eeffdaeff  '  van:J-sW-n€ 

of  the  first  variation  with  respect  to  a'  ,  o'  ,  i '  ,  M  ,  M  ,  M  , 

x  y  xy  x  y  xy 

u  ,  v  ,  and  w  is  then  written  as 


6U"  =  t 


f  J  26°x  (u,x  +  \  V%)  +  2ax  (8u,x  +  W,x5w,x) 


25o '  (  v  +  ^  w2  ]  +  2a  '  f  8v  +  w  6w  | 

y  \  >y  2  ,y  /  y  \  ,y  ,y  ,y  J 


(Continued) 


+  25t  '  (u  +  v  +ww  )  +  2t  '  (8u  +  5v 

xy  ,y  ,x  ,x  ,y'  xy  x  ,y  ,x 


5M  ^ 

+  8v  w  +  w  Bw  ) - w  -  —  8w 

,x  ,y  ,x  ,y'  t  ,xx  t  ,3 


5M  M  8M  M 

— v  -  Sv  +2  — w  +2  -=2-  8w 
t  >yy  t  ,yy  t  >xy  t  ,xy 

f  r  f  f 


'cF'  SF'  cF' 

-  So '  +  - So '  +  -  6t  ' 

So'  x  So'  y  St'  X1 

v  v  wr 


5M  +  -  8M 

SM  y 
y 


+  - SM  dxdy  =  0  (3: 

SM  yJ 

xy 

With  the  grouping  of  like  terms  and  after  integration  by  parts.  Equa¬ 
tion  (32)  becomes 

¥2  f  [/  i  a 

SU"  =  2t  /  /  u  +  -v* - 6o' 

i/a  i/a  l'**  2  ’ 

/  1  2  *'\  /  SF'  \  , 

+  v  +  —  w  -  -  lSo'+lu  +v  +ww - 1  St 

\'y  2  ,y  y  V'y  'x  >x’y  *'/  ** 


w  (S’ ' 
-jt£L+  - 

atf 


/w  dF  \  /  w  dF  \  /w  cF  \ 


w  cF ' 

jxy  _ _ 


-  (o'  +t'  )  8u  -  (o'  +  t'  )  8v 

v  x,x  xy,y'  y,y  xy,x' 


[f°xw,x^x  +  ^°yW^y^y  +  ^TxyW  ,y^x  +  ^Txyw  ,x\ 
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1 

-  (M 

2t  _ 


x,xx  +  My,yy  "  2Mxy,xy)]  Bw 


dxdy 


+  2tf  / 

x  =  +  L/2 

cJSu  | 

f 

dy+2t  / 

t'  6u 

y  =  +  b/2 

-bJ/2 

x  =  -  L/2 

1 J 
-L/2 

xy 

'C 

ii 

i 

o' 

l\3“ 

dx 


+  2t, 


T 

y  =  +  b/2 

c'8v  | 

f 

dx  +  2t„  / 

T  '  &V 

ii 

>< 

/2 

y  =  -  b/2 

-L/2 

xy 

i 

ii 

dy 


X 

X 


l2tf^,x  +  a»V,7  +  V  -  vls>  I  *  ‘  *  ^ 


[2tf^,v  +  *  V  -  V>  I  l  ][  *'l 


dy 


dx 


b/2 


-b/2 


x  =  +  L/2 

Mfiw  | 

,x  x  =  -  L/2 


dy  - 


L/2 


w 

-L/2 


M  8w 
xy  ,x 


y  =  +  b/2 
y  =  -  b/2 


/ 

-L/2 

b/2 


M  5w 

y  ,y 


dx  + 


-b/2 


M  Sw 
xy  ,y 


y  =  +  b/2 
y  =  -  b/2 

x  =  +  L/2 
X  =  -  L/2) 


dx 


dy 


0 


(33) 


For  Equation  (33  )  to  vanish  for  arbitrary  states  of  stress  and  dis¬ 
placement  consistent  with  the  prescribed  displacement  boundary  con¬ 
ditions,  then  each  of  the  terms  must  vanish  identically.  The  Euler 


f 
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equations  and  boundary  conditions  that  result  are  obtained  as  follows: 


otress-displacement  relations: 


do' 

y 


u  + 

,y 


V 


w  w 


Moment-curvature  relations: 


>  (34a) 


'  l 

dM  2t_  ,XX 

x  f 

ap'  1 

2t_  ,yy 

y  t 


ZF'  1 


In-plane  equilibrium  equations: 


)  (34b) 


J 


a'  +  t'  =  0 
x,x  xy,y 


(35a) 


a'  +  t '  =  0 

y,y  xy,x 


(35b) 
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Out-of-plane  equilibrium  equation: 


2t-  [  a  "w  +  2t  '  w  + 
f  x  ,xx  xy  ,xy 


a  w  ] 

y  ,yy 


+  M 


x.xx 


2M  +  M 

xy,xy  y>yy 


0  (36) 


Stress  (natural)  boundary  conditions: 


t'  (x,  ±  b/2)  =  0 


corresponding  to  arbitrary  u-displacement 
along  edges  y  =  ±  b/2 


r'Ji  L/2,y)  =  0 


corresponding  to  arbitrary  v-dis placement 
along  edges  x  =  ±  L/2 


\(±  L/2,y)  =  0 


corresponding  to  arbitrary  rotations  w 
along  edges  x  =  ±  L/2 


My(x,  ±  b/2)  -  0 


L/2 


-L/2 


Oy(x,  ±  b/2)  dx  =  0 


corresponding  to  arbitrary  rotations  w 
along  edges  y  =  ±  b/2 

corresponding  to  constant  v-displacement 
along  edges  y  =  ±  b/2 

(37) 


The  Euler  equations  and  boundary  conditions  established  from  the  vari¬ 
ational  Equation  (33)  are  the  same  as  those  obtained  in  Reference  1  from 
potential  energy  considerations,  except  for  the  addition  of  the  stress- 
displacement  and  moment -curvature  relations  given  in  Equations  (34). 
Naturally,  when  the  constitutive  law  for  the  secant -modulus  deformation 
theory  of  plasticity  (see  Reference  l) 
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is  introduced  into  the  Reissner  functional  [Equation  (26)]  prior  to 
carrying  out  the  variation,  the  functional  can  be  reduced  to  the  po¬ 
tential  energy  function  of  Reference  1.  The  results  obtained  from  the 
variation  are  then  identical  to  those  presented  in  Reference  1. 


APPENDIX  II 

MODIFIED  REISSNER  FUNCTIONAL  FOR  EIASTIC  PROBLEM 


The  modified  Reissner  functional  given  by  Equation  (18),  after  substi¬ 
tution  for  the  assumed  stresses  [Equation  (22  )]  and  displacements 
[Equation  (23)]  and  subsequent  integration  over  the  plate  area,  becomes 


U' 


EV 


=  -  e 


Ap  A~ 

An  +-=■+  — 

3  5. 


+  k„ 


1-2  -2  -2  \ 
(vll+  ^31  +  w13 ) 


4* 


gl  +  S=fi)  3  .  . 


-2 

V13 

^  ull”l3  +  k34  -f 


+  A- 


(7U+^l) 


+  k35wnwl3  +  k3g 


-2 

W13 


+  C 


^1  -2 

M*5  T  +  k6w3l  -  Vnv31  )+  vllvi3k37 


S2 

+  ,  k38  -  "ll”3A5 


+  k 


■8 


“2  n 

-2  -2  W3l  ^  — 

W11  +  ^13  +  “  +  “  vllv3l 


+  Bc 


(vll  +  9^13  +  ^uv31) 


+  klowllw3l  +  kll 


-2 

w3l 


+  c 


ki2  /  -2  w3i  3  _  _  \ 

n2  ^  wn  +  —  +  -  W11W31  j  +  k39wl3 


(Continued) 
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+  Vn  +  V'n  +  (39> 

where 

7n  =  (r)vn 
w3i  =  (r)v3i 
W13  =  (f")w13 

The  coefficients  k1  through  k^  are  numerical  constants  resulting 
from  the  integration  process.  The  constants,  some  of  which  depend  upon 
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either  or  both  of  the  parameters  b/X  and  ith/b  ,  are 


kl 

= 

2.46739  [p2/(jth/b f] 

kl7 

= 

0.01975 

k2 

= 

0.13069 

kl8 

= 

0.08889 

k3 

= 

0.04110 

kl9 

= 

0.07111 

h 

= 

0.06219 

k20 

= 

0.00054  6* 

b 

= 

0.35653 

k21 

= 

0.15238 

k6 

= 

1.05684 

^2 

= 

0.00931  p2 

*7 

3 

0.39268 

3 

0.00798  P2 

k8 

3 

2.46739  tl/(jth/b  f] 

*24 

3 

0.00117  P2 

k9 

3 

-0.15890 

^5 

= 

0.41123  Pk  +  2p2  +  1 

kio 

= 

0.08592 

k26 

= 

0.41123  (8lpk  +  1802  +  1) 

kn 

= 

-0.04199 

^7 

= 

0.00035  P2 

k12 

= 

0.09019  p2 

k34 

= 

0.31082 

kl3 

= 

0.31831  [p2/(rth/b f] 

k35 

3 

-0.07298 

kl4 

= 

0.05939 

k36 

= 

0.15801 

k15 

= 

0.00276 

k37 

= 

-0.0108 

kl6 

= 

0.17780 

k38 

= 

O.OI692 

(Continued) 
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■39  =  0.09069  6" 

k43  =  0.19^90 

;4l  =  -0.00526  B2 

k44  =  0.^1123  (P4  +  18  S2  +  81) 

:4l  =  0.04533 

k45  =  -0.01192 

%2  =  -0.05797 

(41) 

The  requirement  of  the  vanishing  of  the  first  variation  of  the  functional, 
given  by  Equation  (39)  with  respect  to  the  arbitrary  parameters,  leads  to 
a  set  of  nonlinear  algebraic  equations,  the  solution  of  which  establishes 
the  values  of  the  stress  coefficients  A^  ,  Ag  ,  A^  ,  Bg  ,  B^  ,  and  C 
and  the  displacement  coefficients  w^  ,  w^  ,  and  w^  .  The  results 
obtained  with  =  0  are  designated  Solutions  Al,  whereas  the  results 
obtained  with  w^  retained  are  designated  Solutions  A2. 


APPENDIX  III 

UNMODIFIED  RE  ISSUER  FUNCTIONAL  FOR  EIASTIC  PROBLEM 


The  unmodified  Reissner  functional  for  elastic  problems  given  by 
Equation  (1*0,  after  substitution  for  the  assumed  stresses  [Equations 
(22)  and  ( 24) ]  and  displacements  [Equation  (23)]  and  subsequent  in¬ 
tegration  over  the  plate  area,  becomes 


U" 

EV 


-  e 


Ap  Ao 

At  +-+-*- 

1  3  5 


+  k. 


-2  2  —2 
w  +  9V-3J+W 


+  A, 


3,  ,  4 


-2 

13 


*2  “T  ^  “  72  vnvi3  +  “f  k34 


+  A. 


2jt 


31 1  _  _  wi3 

2  +  k35Wllv13  +  ~  k36 


(=11  +  9&) 


+  c 


+  k6=3i  *  Vn=3i)  +  Wu 


"l 

+  k38  Y  -  Vl3"31 


+  k 


8 


B_ 


-2  -2  W3l  ^  — 

W11  +  ^13  +  9  +  2  WHV 31 


+  B, 


(^1+  & 


k  ;  ;+  k  3i 

2  *10  11  +  *11  2 


(Continued) 
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“ 

k12| 

-£  W31  3  -  -  ^ 

l  .  1  .  ,  ( 

2 

n 

^11  9  +2WllW3ly 

r  k39v13  +  Ho  l 

+  c 


4.  —  v  W 

+  2  13  31 


-  Ck 


13 


( -e  w3i  ^  —  \  -e  _  _ 

kl4  [vn  +  -f  +  -  W11V31  j+  k4iwl3  +  k42wnwi3 


+  k43wl3v3l 


-  -  p  An  (wn  +  a  wi3)  -  -  bu  (vu  +  9  «  wl3 ) 


+  -  cn  P  (*u  +  3  a  vl3) 


1 

2 


A2t^+^tC2(ki5+ll£0+^-,^7) 


(^) 


(a  11  +  Bn  ”  AnBu  +  3cn^ 


222 

+  B2kl8  +  B3k19  +  -  AjAg  +  -  AXA3  +  -  A^ 


-  C  (Agkl6  +  A3k17  +  Bgkgg  +  B3kg3)  +  BgB^j^ 


(42) 


where  ,  w31  ,  and  vl3  are  given  by  Equations  (4  0)  and  the  k  co¬ 
efficients,  by  Equation  (4l). 
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The  results  obtained  from  the  simultaneous  vanishing  of  the  first  vari¬ 
ation  of  the  functional  given  in  Equation  (42)  with  respect  to  the 
stress  coefficients  A1  ,  Ag  ,  ,  B^  ,  ,  C  ,  A^  ,  ,  and 

a  and  the  displacement  coefficients  ,  and  are 

designated  Solutions  B. 


APPENDIX  IV 

NEWTON -RA PHS ON  ITERATIVE  TECHNIQUE 


The  requirement  that  the  Reissner  functionals  given  in  Equations  (39), 
(i^2),  and  (25)  have  a  stationary  value  results  in  sets  of  simultaneous 
nonlinear  algebraic  equations,  the  solution  of  which  establishes  the 
values  of  the  stress  and  displacement  coefficients.  The  basic  Newton- 
Raphson  technique  is  employed  to  effect  the  solution.  The  method  is 
illustrated  by  the  following  example.  Consider  a  function  F(x^,Xg, 

•  •  •  •  xffl)  for  which  the  values  of  x^  are  sought  that  render  F 
stationary.  This  requires  that 


aF(x±) 

-  =  ♦i(xi)  =  0  i  =  1,2,  ...m  (number  of  independent 

variables ) 


If  approximate  values  of  are  known,  the  above  relations  may  not  be 
satisfied.  Therefore,  by  applying  a  suitable  correction  A  x^  such  that 
♦^(x^  +Axi)  *  o  ,  the  desired  values  of  x ^  are  obtained.  To  deter- 
mine  the  correction  A  x  ^  ,  the  functions  ^  (x^  +  ax^)  are  expanded 

into  a  Taylor  series  about  xJ  with  only  the  linear  terms  retained. 

This  procedure  results  in  the  following  system  of  linear  algebraic 
equations 


♦l(*i 


3 


o 


(U3) 


where  i  =  l,2,3...m  .  Written  in  matrix  form,  the  equations  are 


d2 F 

fc1  , 

Ax. 

■  - 

*1*5. 

i 

*i 

(WO 


where  1  =  1,2... m  and  J  =  l,2...m  .  Thus,  by  adding  the  correction 

Ax^  to  the  starting  value  Xj  and  by  repeating  the  process  until 
cF/cb^  -*  0  ,  the  solution  is  obtained. 


For  the  inelastic  problem  formulated  in  Equation  (25),  the  partial 
derivatives  indicated  previously  are  obtained  numerically,  whereas, 
in  the  elastic  solutions,  the  derivatives  are  evaluated  explicitly. 
The  partial  derivative  approximations  utilized  are 


cJF  F^  +  &xA)  -  F^  - 

■  —  S  I  «  ■  — 

cfcj  25x1 


b^F  F(xa  +  6xi)  -  2?^)  +  F(xi  -  bx±) 

dx*  (S^)2 


b^F  F(x1  +  5x1  ,  Xj  4-  &Xj  )  -  F(x,^  +  bx^)  -  F(x^  4-  6x^  )  +  F(x1,Xj) 

(^5) 

where  i  ■  1,2... m  and  J  =  l,2...m  .  - -  ■  ■■ 

The  step  sizes  6x.  used  in  evaluating  these  derivatives  are  taken  as 
-3  1 

x±  x  10  .  The  area  integration  indicated  in  Equation  (25)  is  per¬ 

formed  numerically  utilizing  a  basic  trapezoidal  technique  carried  out 
on  a  5  x  5  plate  grid. 
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